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the anticrossings for resonances involving the lower two levels in regions away from 
accidental or low-order resonances of the upper two levels. 
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1. Introduction 

In early studies of spin dynamics in magnetic fields, Bloch and Siegert considered 
the basic problem of the two-level system with a sinusoidal perturbation pQ. The 
perturbation increases the level separation (the Bloch-Siegert shift), and energy can 
be exchanged with the two-level system when the transition energy matches an odd 
multiple of hu>o (Bloch-Siegert resonances), where uo is the oscillatory frequency. These 
effects have been studied with a dynamical Hamiltonian [HE] (the Rabi Hamiltonian) , 
and also with a static Hamiltonian |3j (the spin-boson Hamiltonian) in which the two- 
level system and oscillator are modeled as coupled quantum systems. These models 
have been of interest over the years for applications to physical systems, including spin 
problems [H EU [5] and atoms in strong electromagnetic fields [21 El IE] ; and also as 
standard model problems on which new approximation methods can be tested [HI [10] . 
We have been interested in such problems in order to better understand energy exchange 
between quantum systems with highly mismatched characteristic energies. 

In this work, we address a generalization of the spin-boson problem in which a three- 
level system is coupled to an oscillator. Our interest in this problem is focused primarily 
on the multiphoton regime, in which the characteristic energy of the oscillator is small 
compared to the available transition energies of the three-level system. While there 
exists a modest literature on coupled oscillator and three-level models [TT 1 [T2 l [T3l ITU 115] . 
there does not appear to be available previous work relevant to the multiphoton 
region. The rotating wave approximation has been applied with success near low-order 
resonances [HI [TTJ, [TEJ HH EQl EH 122], but one would not expect this approach to be 
helpful in the multiphoton region. 

The replacement of the two-level system of the spin-boson problem by a three- 
level system in the generalization of the model considered here very much complicates 
the problem. The convenient mathematical machinery of the spin operators and Pauli 
matrices available for the spin-boson problem is now replaced by the less convenient 
SU(3) operators and Gell-Mann matrices [23], and the diagonalization of the three-level 
system now involves a cubic characteristic equation. There can occur in addition level 
crossing effects in the three-level system in which the middle level is pushed through 
one of the other levels. Faced with such difficulties, it is no surprise that the three-level 
version of the model has not received comparable attention as a workhorse problem in 
comparison with the simpler spin-boson problem. 

Here, we propose to use the more complicated coupled oscillator and three-level 
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model to test an approximation scheme that appears to have worked well in the spin- 
boson problem [24J, and also in the generalization to the spin-1 version of the spin-boson 
model [25]. In this approach, the two- level or three-level system is diagonalized directly, 
treating the oscillator operators parametrically. The rotation which accomplishes this 
produces weaker residual couplings when applied to the oscillator. This rotation was 
discussed early on for the spin-boson problem by Wagner [26J. In the multi-photon 
region, the energy levels of the full problem are well approximated away from resonance 
by simple averages over the parameterized diagonalized energy levels. We expect this 
also to be the case in the coupled oscillator and three-level system problem. In the 
first half of this paper, we develop the approximation explicitly and compare with exact 
results for a test problem, with the result that the approximation gives good results away 
from resonances. We found previously in the spin-boson model that we could understand 
the level splittings at resonance from a simple two basis state approximation using matrix 
elements of the residual interaction arising from the rotation of the oscillator. We are 
interested in whether this is also the case in the more complicated version of the problem 
under consideration here. In the second half of this paper, we develop the interaction 
terms explicitly and compare approximate level splittings with exact numerical results 
from a test problem. We find that the method works well for high-order resonances 
away from accidental or low-order resonances on other transitions. 
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Figure 1. Three-level system and coupling. 



2. Basic model and unitary transformation 



The coupled oscillator and three-level system, which generalizes the spin-boson model, 
can be described using the Hamiltonian 
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In this model, the three-level system has unperturbed energies E X) E 2 , and E 3 . 
Transitions between the first two levels are described by a linear coupling term, where 
the interaction has strength U . Transitions between the second and third levels are 
also included through another linear coupling term, with interaction strength V. Such 
a coupling scheme is normally called a ladder or a H configuration [27]. We will also 
assume that 



\E{ — Ej\ ^> Kujq , n 3> 1 



2.1. Rotation 

In our earlier work, we used a unitary transformation to rotate the spin-boson 
Hamiltonian in to a more complicated form which we wrote as 

iftHU = H + V + W (2) 
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The rotation diagonalized the two-level model parameterized by the oscillator variable 
y defined as 

a + a) 



y 



V2 



(3) 



Although the transformation introduced terms more complicated mathematically, the 
rotated Hamiltonian was conceptually simple. The first term H was found to produce 
a reasonably good approximation to the energy levels away from the anticrossings. The 
second term V can be used to estimate the level splittings at the anticrossings. The 
third term W is small in the large n regime, as is the case here. 



2.2. Unperturbed rotated Hamiltonian H 

We focus first on H with the goal of developing an approximation for the energy levels 
away from the anticrossings. Following the approach used in the spin-boson version of 
the problem, the unitary transformation we use here is one which diagonalizes the part 
of the Hamiltonian that does not include the oscillator 

H -hu ata = U^[h - Uuqo) a^jU (4) 

Combining the different 3x3 matrices that appear in the Hamiltonian, we may write 
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2.3. Diagonalization 



An explicit construction of the unitary transformation could be done for the simpler 
spin-boson problem, but this is inconvenient for the three-level version of the problem 
that we are interested in here. Nevertheless, we can still make progress by noting that 
the rotation accomplishes a diagonalization, which can be done by solving 
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Figure 2. Energy levels Ej(y) (in units of Hujq) as a function of y for a model with 
Ei = 0, E2 = 11 hujo and £3 = 24 huio. The number of oscillator quanta n is 10 8 . (a) 
The coupling strengths are U^/n — 0.8 and Vy/n = 0.8. (b) The coupling strengths 
are Uy/n = 0.1 and V yfn = 1-0. 



This leads to a characteristic equation that is cubic 

(E 1 - E)(E 2 - E)(E 3 -E)- (E 1 - E)2V 2 y 2 - (E 3 - E)2U 2 y 2 = (7) 

We denote the solutions to this cubic equation as Ei(y), E 2 (y) and E 3 (y), defined so 
that 



E 1 (0)=E 1 



E 2 (0) = E 2 



£3(0) = E 3 



The solution of the cubic equation and expressions for the eigenvalues Ej(y) are discussed 
in Appendix A. The unperturbed part of the rotated Hamiltonian Hq is then 

I E x {y) \ 
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2-4- Level crossing and anticrossing 

The diagonalization of the three-level system parameterized by y produces energy levels 
that depend on y, examples of which are illustrated in Figure El We expect static 
coupling between levels to push them apart, and in the calculation of Figure [2](a) this 
is exactly what we see. However, since the interaction pushes levels apart, it is possible 
that one level may be pushed into and through another level. This is illustrated in 
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Figure 12(b), where the interaction between levels 2 and 3 is strong, and level one 
weakly coupled. In this case level 2 anticrosses level 1. We refer to the middle level 
E 2 (y) on either side of the anticrossing. 
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3. Energy levels away from anticrossings 

To understand the systematics of the energy levels we consider solutions to the time- 
independent Schrodinger equation 

= (9) 

Some of the energy levels from such a computation are illustrated in Figure EJ where we 
have defined the dimensionless coupling strengths g\ and g 2 according to 

91 = p 22 = p r- (10) 

There are many energy levels; we expect three levels for each n, which is the same as 
we would get in the absence of coupling. As one can see from Figure [3J the complete 
spectrum is locally made up of three curves with the same pattern repeated over and 
over again with increasing energy (where the energy offset is 2 hu>o). The actual offset 
is huo, but we have only included the even parity (even total index j plus n) energy 
levels to make the figure simpler. The energy levels of the odd parity (odd total index j 
plus n) are nearly identical, shifted by one unit of hu . This periodicity occurs because 
in the large n limit, the coupling of n + 1 quanta is little different from coupling of n 
quanta. Under these conditions, the energy levels can be parameterized according to 

E j , n (g u g 2 ) = E j (g 1 , g 2 ) + nhu (11) 
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The corresponding physical statement is that the oscillator impacts the three-level 
system strongly, causing the levels to shift; while the three- level system has very little 
impact on average on the oscillator, so that its levels are not shifted. This situation is 
similar to that observed in the spin-boson problem when n is large [23] . This motivates 
us to try to understand the shifted (or dressed) energy levels Ej(gi,g 2 ), as these are 
fundamental to the model in this limit. 

3.1. Dressed energy levels from the rotated Hamiltonian Hq 

We can develop an effective estimate for the energy levels away from the anticrossings 
using an approach similar to that used for the spin-boson problem. We base our 
approximation on the eigenvalues of the rotated H problem which can be written as 

/ E x {y) \ 
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We can separate the three-level degrees of freedom from the oscillator degree of freedom 
by using wavefunctions of the form 

(1\ (<>\ (*\ 
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The Uj(y) functions satisfy a one-dimensional Schrodinger equation of the form 



Mv) (is) 



E + ^fu^oj Uj(y) = Ej(y) + ^hu 



(14) 



This is similar to what we obtained in the spin-boson problem [24] . 
3.2. Variational estimate 

A useful approximation valued away from the level anticrossings can be developed by 
adopting trial wavefunctions in which the Uj(y) function is taken to be a simple harmonic 
oscillator functions 



My) = ^niv) 



(15) 



Such a wavefunction is relevant in the large n limit for a variational estimate of the 
energy, leading to the approximation 



Ej,n{gi,92) = (n\Ej(y)\n) + nhu = Ej(gi_, g 2 ) + nhu Q 



(16) 
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Away from the anticrossings, the eigenf unctions of the original Hamiltonian H have 
average n and j values which are close to the integer values which we assign in the rotated 
Hq problem here. When n is large, we can take advantage of the WKB approximation 
to write 

(n\E 3 {y)\n) = - f -^M=dy (17) 

71 J-e \Jt — y 2 

with e = 2n + 1. When n is large, the WKB approximation gives excellent agreement 
with the results from a brute force solution of Equation (TH1) . as was the case previously 
for the spin-boson problem. An algebraic expression in the case of the middle level is 
given in the Appendix. 
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4. Model test problem 

To illustrate the approximation under discussion, we consider a specific example. In 
this example, we are interested in the basic question of whether the approximation 
described in the previous section is a good one away from resonance. However, since 
the three- level model and coupled oscillator problem is considerably more complicated 
than the analogous spin-boson problem, our task becomes more difficult. We focus on 
a three-level system with unperturbed energy splittings given by 

E 2 -E l = 11 huo E 3 -E 2 = 13 huj (18) 

with gi and ranging between and 1, and g 2 ranging between and 1.2 (we extend the 
range for g 2 since interesting features extend out a bit further along the g 2 axis). 

4-1. Resonances associated with transitions between levels 1 and 2 

We expect that the approximation for the energy levels discussed above will be best 
away from resonances, so our first goal is to understand where resonances occur. The 
most important resonances are those involving transitions between levels 1 and 2, and 
levels 2 and 3. In the case of resonances involving the lowest two levels, the resonance 
conditions is given approximately by 

E 2 (gi,g2) - E 1 (jg 1 ,g 2 ) = Anhu (19) 

where An is odd. The associated contours are illustrated as a function of the 
dimensionless coupling constants g\ and g 2 in Figure HJ based on the WKB 
approximation. When both g\ and g 2 are small, it should take 11 oscillator quanta 
to match the transition energy since the energy levels in this case are only weakly 
perturbed. For small g 2 , the dressed transition energy E 2 (g 1 ,g 2 ) — Ei{gi,g 2 ) increases 
with increasing g 1; so that more oscillator quanta are required for a resonance. This can 
be seen in Figure H] in the resonance contours for increasing An coming down to the g\ 
axis. In the upper left of this figure, things are much more complicated. As g 2 increases, 
the splitting between the upper two levels is increased. If g\ is small, then level 1 is only 
weakly coupled, and it is possible for level 2 to be pushed near level 1. The separation 
between the lower two levels is decreased, so that lower-order resonances (in which An 
is smaller) occur. The coupling due to the V operator between the two lowest levels 
is much stronger for these lower-order resonances, so that the mixing is much greater 
and the associated level splitting is much larger. Approximating these states as pure 
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Figure 4. Contours for resonances between levels 2 and 1 (from the rotated Hq 
Hamiltonian), with the exchange of an odd number of oscillator quanta. The number 
of oscillator quanta An is indicated with each resonance. 
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Figure 5. Resonances on transitions between levels 1 and 2 through the inverse 
magnitude of the resonance condition |£2(<7i, 32) — E\(g\, 52) — AnUuJo^ 1 as a function 
of gi and g 2 ■ Computations for this plot are done using the original unrotated 
Hamiltonian H. The results are in units of (Hujo)^ 1 . 



eigenfunctions of the Hq Hamiltonian becomes a much poorer approximation in this 
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regime. 

We can see the resonances in the results shown in Figure [5] taken from a direct 
computation using the original H Hamiltonian. One can see that most of the resonances 
predicted by the WKB approximation appear in the full computation. 

4-2. Resonances associated with transitions between levels 2 and 3 

The resonance condition associated with transitions between levels 2 and 3 is given 
approximately by 

E 3 (gi, g 2 ) - E 2 {gi, g 2 ) = Ank (20) 

where once again the number of oscillator quanta exchanged must be odd. This 
resonance condition from WKB calculations is illustrated in Figure [61 For small g\ 
and g 2 it now takes 13 oscillator quanta to match the transition energy, since in this 
case it is only weakly perturbed. When g\ is small, an increase in g 2 pushes levels 2 
and 3 apart, increasing the number of quanta required for a resonance. This situation 
is qualitatively similar as considered above, except that the axes are reversed. When 
g 2 is small, large values of g\ can cause levels 1 and 2 to push apart strongly, driving 
level 2 toward level 3. Due to the asymmetry between the level splittings, this occurs 
for gi values larger than those for g 2 in the previous case for resonances between levels 
1 and 2. The resulting small separation that occurs in the lower right hand part of 
Figure [6] between levels 2 and 3 leads to the occurrence of low-order resonances, with 
the associated strong mixing and large level shifts due to the V operator. 

Once again we show analogous results from direct calculations based on the original 
H Hamiltonian in Figure [71 Most of the resonance lines from the WKB calculation are 
evident in the full H calculation. The effects of the low-order resonances between levels 
3 and 2 are apparent in the lower right part of the plot. 

4-3. Energy levels 

From the discussion above, we expect that the WKB energy level estimates based on 
the rotated Hamiltonian H should be good away from regions where resonances occur, 
and should be best where the separation between the nearest level is largest. In Figure El 
we illustrate the magnitude of the difference between the exact energies level 1 obtained 
from the original H problem, and the WKB approximation: 

«,<»,«,) -i/* ^MLiy 

TT J-e We — y z 
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Figure 6. Contours for resonances between levels 3 and 2 (from the rotated Hq 
Hamiltonian), with the exchange of an odd number An of oscillator quanta. The 
contour labels indicate An. 
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Figure 7. Resonances on transitions between levels 1 and 2 through the inverse 
magnitude of the resonance condition |-E 3 (<7i, g 2 ) — E 2 (gi, 32) — Anftwol -1 as a function 
of gi and g 2 ■ Computations for this plot are done using the original unrotated 
Hamiltonian H. The results are in units of (Hujo)^ 1 - 



One sees that when gi and #2 are small, that the difference is also small, and the 
approximation is very good. In the upper left part of the plot there occur the largest 
deviations as we expected due to low-order mixing effects not included in the Hq 
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Figure 8. Magnitude of difference between the exact result for i?i (51, 52) 

from numerical calculations of the unperturbed Hamiltonian H, and the WKB 
approximation for E\(g\,g2) in units of Hloq. 

problem. Elsewhere the agreement is generally rather good away from the resonances. 
The energy shift associated with the high-order resonances is not so large, and deviations 
are observable only when gi becomes large. We note that the energy level separation 
between the lowest two levels away from the upper left part of the plot is between 
about 10 and 40 hu ; hence the deviation between exact and approximate results is 
generally less than 0.1% of this separation. In the region where we would expect the 
approximation to be relevant, we find that it is very good. 

In Figure [9] we show similar results for the magnitude of the difference in energy 
for level 3. Once again, the difference is small when gi and #2 are small. Away from the 
lower right region (where level 2 is pushed toward level 3) the differences are small and 
the approximation is good. In the lower right region, strong mixing occurs due to the 
low-order resonances, and the approximation under discussion degrades. 

In the case of level 2, the region where the approximation is accurate becomes 
further restricted, as can be seen in Figure [10J The strong mixing that we previously 
encountered in the upper left part of the plot in the case of level 1 is reproduced in the 
case of level 2, since level 1 is mixing primarily with level 2. The strong mixing associated 
with level 3 in the lower right part of the plot in the case of level 3 is also reproduced 
in the case of level 2, since level 3 is mixing with level 2 as well. Consequently, the best 
results are obtained in the central region, away from the problem regions in the upper 
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Figure 9. Magnitude of difference between the exact result for E% (31,92) 
from numerical calculations of the unperturbed Hamiltonian H, and the WKB 
approximation for E3 (gi, 92) in units of Hu>q. 
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Figure 10. Magnitude of difference between the exact result for £2(91,92) 
from numerical calculations of the unperturbed Hamiltonian H 7 and the WKB 
approximation for £2(31,92) in units of Tiloq. 



left and lower right part of the region. 
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5. Remaining terms of the rotated Hamiltonian 



We set out initially to implement a unitary transformation similar to the one which we 
studied previously in the spin-boson problem, in order to obtain a rotated Hamiltonian 
of the form 



hPEU = H + V + W 



In previous sections we developed the first term in the rotated Hamiltonian H , and 
compared results from computations with those derived from the original Hamiltonian 
H. To complete the unitary transformation, we need to rotate the simple harmonic 
oscillator part of the Hamiltonian to obtain V and W . We may write 



V + W = hP 



huJocPa 



(21) 



5.1. Transformation under rotation 

Perhaps the simplest way to approach this problem is to examine how the simple 

harmonic oscillator operators transform under the unitary transform. The oscillator 

operators under discussion here are y, and — , where 

dy 



d_ 

dy 



a — d^ 



Since the rotation diagonalizes the three-level system for each y, the unitary operators 
itself depends on y. We may write 



U = U{y) 

to exhibit this explicitly. There is no dependence of the unitary operator on 
Consequently, the rotation of the oscillator operators can be written as 

U\y)yU{y) = y 



y 



dy 



d_ 

dy 

(22) 
(23) 



5.2. Rotation of the harmonic oscillator 



We can make use of this transformation to construct the other terms in the rotated 
Hamiltonian. We have 
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Inserting expressions for the transformed operators produces 



a) a 
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This allows us to identify the two remaining parts of the rotated Hamiltonian, which 
we may write as 



V = 
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(26) 
(27) 



5.3. Reduction ofU ] {y)^-U{y) 

It is possible to simplify these operators by taking advantage of properties of the unitary 
matrix. The unitary operator satisfies 



U\y)U{y) = 1 

independent of the parameter y. Hence, if we differentiate in y we obtain 
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dy 
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(29) 



A consequence of this is that the matrix which represents W(y)—U(y) is antisymmetric 
if we chose U (y) to be a real matrix. In this case we may write 



/ U u U 21 Z4i \ 

IA\2 U-22 1^32 
\ ^13 ^23 ^33 / 



( U' u U[ 2 U{ 3 \ 

^21 ^22 ^23 
\ ^31 ^32 ^33 / 



/ 



F 13 \ 



F12 ri 3 
F\2 F23 

\ ~ i*13 -^23 J 



(30) 



where 
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The Fij matrix elements can be obtained directly from the unitary matrix, which allows 
us to write 



F±2 


= UwU'yi + IA21W22 + 


(31) 
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The U\y)^j-U{y) operator can then be expanded as a sum of spatial parts and SU(3) 
matrices according to 
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5.4- Reduction ofV 

We can use these results to recast the V in a similar form. We may write 



V 



( -i \ 

i 
v j 
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(35) 





\ i y 

/ \ ' 

7 - 

+ ^F 23 {y) + F 23 (y)— -i > 

The y operator written in this form is seen to be closely related to the one obtained in 
the spin-boson problem. A similar reduction of W can be carried out without difficulty 
from these results. 
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6. Level splitting at the anticrossings 

As discussed above, the energy levels away from the anticrossings are given 
approximately by 

E j>n = Ej{g x ,g 2 ) + nhu 

The approximate resonance conditions for the level anticrossings between the lowest two 
level as discussed above is 

E 2 (gi,g 2 ) - E 1 (g 1 ,g 2 ) = Anhu 

for An odd. When we studied the rotated spin-boson system using a similar unitary 
transformation, we found that the level splitting was well approximated using degenerate 
perturbation theory in the rotated frame based on matrix elements of the V operator. 
We expect a similar behavior in the three-level version of the problem, which would 
result in an estimate for the level splitting given by 

5E = 2|($ 1>n |V r |$ aift _ An )| (36) 
6.1. Level splittings, large An 

As discussed in Section 4, the regions with different g% and g 2 would be expected to 
behave differently due to the presence of low-order resonances. Hence, we begin the 
discussion with a comparison in a relatively benign region in which levels 1 and 2 have 
large separation, and we stay away from the problem area in the lower right part of the 
plots discussed in Section 4. We seek resonances along a line with 

g 2 = 0.3 gi 

Along this line the resonances between the lower two levels do not overlap other 
resonances; the resonances are all of high order, and there is no interference from the 
problem area associated with levels 2 and 3. Consequently we would expect the lower 
two levels to give good agreement, much as we found in the spin-boson problem with a 
two-level system. Results from computations based on the model test problem discussed 
in Section 4 are shown in Figure [TT1 We see in this figure very good agreement, similar 
to what we found in the spin-boson case |24j. 
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Figure 11. Level splitting at the anticrossings between levels 1 and 2 for the model 
described in Section 4, with gi = 0.3g±. Black circles: degenerate perturbation 
theory based on wavefunctions in the rotated frame; gray circles: results from a direct 
computation using the original unrotated Hamiltonian H. Energies are given in units 
of Tiujq. 

6.2. Interaction with low-order resonances 

We would expect a degradation of the approximation in the vicinity of the region where 
the other transition experiences low-order resonances. We continue to focus on the 
resonances involving the lowest two levels, this time on a line given by 

92 = 0.1 g x 

which goes into the problem area with low-order transitions between the upper two 
levels. Results are presented in Figure [12j We see in this case good agreement between 
the results from degenerate perturbation theory based on the rotated problem, and 
results from the original H problem everywhere except at resonances with An equal 
to 23 and 25. At these points the upper two levels are mixed due to strong low-order 
V interactions, and one finds two resonances in the vicinity of the single resonance 
predicted by the rotated H model. 

6.3. Other effects 

We find poorer agreement when we move into the region in which g 2 is greater than g\. 
We consider in Figure [13] resonances along the line defined by 
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Figure 12. Level splitting at the anticrossings between levels 1 and 2 for the model 
described in Section 4, with 52 = O.lgi. Black circles: degenerate perturbation 
theory based on wavefunctions in the rotated frame; gray circles: results from a direct 
computation using the original unrotated Hamiltonian H. Energies are given in units 
of hujo- 



92 = 1.1 9\ 

The results in this case can be understood as being due to two different effects. The 
level splittings at the anticrossings for the full H problem begin to fall systematically 
below those obtained from degenerate perturbation theory for the rotated problem once 
we consider the g% = g\ line. Two of the resonances (13 and 17) are lower still due 
to interference from nearby resonances involving the upper two levels. We find poorer 
agreement between the exact numerical results for the H problem and those from the 
approximation in this region, and worse still as the ratio of g 2 to g 1 increases further. 
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Figure 13. Level splitting at the anticrossings between levels 1 and 2 for the model 
described in Section 4, with gi = l.lgi. Black circles: degenerate perturbation 
theory based on wavefunctions in the rotated frame; gray circles: results from a direct 
computation using the original unrotated Hamiltonian H. Energies are given in units 
of hujo- 
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7. Summary and conclusions 

In our previous work, we considered the unitary transformation of the spin-boson 
problem which resulted in a rotated Hamiltonian that was functionally simpler (although 
more complicated mathematically). One part of the rotated Hamiltonian (H ) seemed 
to give rise to an underlying unperturbed problem in which the energy levels were very 
good, and in which no interaction was present between the different levels. Another 
part (V") seemed to give a good account of the interaction responsible for the level 
splittings at the anticrossings. The remaining part (W) is small in the large n limit, 
and could be neglected. In this work, we considered a generalization of the spin-boson 
model in which a three-level model is coupled to an oscillator. We were interested in 
whether a similar rotation could be developed, and whether it would exhibit similar 
useful properties. It was conjectured in Ref. (28] that the rotation could be applied 
to multilevel generalizations of the spin-boson problem, and the results presented here 
constitute such an example. 

We succeeded in implementing such a rotation, which allowed us to develop 
analogous approximations for energy levels and for level splittings at the anticrossings. 
To study the approach, we developed a model test problem to provide a concrete example 
that could be used to compare results from numerically exact calculations with those 
from the rotated version of the problem. This model problem itself illustrated interesting 
new effects not present in the spin-boson problem, in which strong coupling between 
two levels result in one of the levels being pushed close to the third level. The resulting 
problem areas have no analog in the simpler spin-boson model. 

The energy levels predicted from the rotated Ho problem are in excellent agreement 
with those of the full H model away from the resonances, and away from the problem 
regions in the vicinity of low-order resonances. Approximate energy levels of the 
dressed problem were computed using the WKB approximation; which is convenient 
for numerical calculations, and which is in good agreement with energy levels obtained 
through brute force numerical solution of the H problem. Energy levels were generally 
found to be good to better than 0.1 Kuoq away from the problem regions, which 
corresponds to a relative accuracy of better than 0.1 % compared to the transition 
energies. The accuracy is similar to what is observed in the spin-boson problem for 
similar transition energies and dimensionless coupling strengths. 

The computation of level splittings at the anticrossings requires the construction 
of interaction terms in the dressed problem, which is more difficult than in the spin- 
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boson problem because of complications due to the three-level model. We presented a 
straightforward method which relies on properties of the unitary transformation matrix, 
and can be implemented in practice relatively simply using numerical differentiation. 
Using these results, we found that the level splittings could be determined accurately 
using degenerate perturbation theory (similar to the case in the spin-boson problem); as 
long as the dimensionless coupling strength (g 1 or g 2 ) of the transition in question was 
the larger of the two, and as long as computations are done away from the problem area 
with low-order resonances on the other transition. We expect that improved results 
for the level splittings can be obtained by employing somewhat more sophisticated 
approximations than first-order degenerate perturbation theory. Few-state models in 
the rotated frame that take into account interfering resonances would be expected to 
extend the range over which the level splittings could be approximated accurately. 

The dressed problem resulting from the unitary transformation provides a different 
view of the coupled three-level and oscillator problem which allows us to understand 
the multiphoton regime. The new energy level approximation seems to work well, and 
may be useful for applications. The approximation for calculating level splittings is less 
robust, but also seems to be useful as long as applied in trouble-free regions. We would 
expect this approach to apply to more complicated models involving more levels, which 
would provide similar predictive capability for state energies. In general, the approach 
should be most useful in regimes in which the oscillator energy is small compared to all 
relevant transition energies. 
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Appendix A. Solution of the eigenvalue equation 



The characteristic equation for the energy eigenvalues that we obtained in the 
diagonalization of the three-level system is the cubic equation 

(E 1 - E)(E 2 - E){E 3 -E)- (£1 - E)2V 2 y 2 - (E 3 - E)2U 2 y 2 = 

In this Appendix, we review the solution for the three energy eigenvalues (the solutions 

of which are well known, but not usually presented in the form we discuss). Our strategy 

will be to reduce it to the form of the triple angle sine formula 
6 6 

4sin 3 --3sin- = -sin 6 (A.l) 

o o 

Appendix A.l. Offset 

We first define a new energy variable that is offset by the average energy of the three 
states 

e = E — Ei + ^ 2 + jS3 (A.2) 
which allows us to rewrite the characteristic equation in the form 

e 3 -ae = (3 (A.3) 
where a and (3 are given by 

a = El + El + gj ~ ElE2 ~ ElEs - E2E " + 2U 2 y 2 + 2V 2 y 2 (A.4) 



P = ^(Ef + El + El) - ^[E 2 (E 2 + E 3 ) + E 2 (E X + E 3 ) + E 2 ^ + E 2 )\ 

+ ^E 1 E 2 E 3 + hl 2 y 2 {E l + E 2 - 2E 3 ) + ^V 2 y 2 (E 2 + E 3 - 2E 1 ) (A.5) 

Appendix A.2. Scaling 
Next, we scale according to 

Q 

e = Asm- (A. 6) 

This leads to 

4sm V^ sm 3 = T» (A ' 7) 

The triple angle sine formula is recovered with the identifications 
4a 

# = 3 ( A - 8 ) 
43 

= -sin9 (A.9) 



Multiphoton Block- Siegert shifts and level-splittings in a three-level system 27 
Appendix A. 3. Energy eigenvalues 

Since the sine function is invariant under shifts of multiples of 2n 

sin6 = sin(# + 27r) = sin(# + 4vr) (A. 10) 

there are three solutions to the triple angle sin formula 

Consequently, we obtain three solutions to the eigenvalue equation which we may write 

as 

Ei + E 2 + E 3 . / 9 \ . . 

£ = + Asml-j (A.ll) 

Ei + E 2 + E 3 . /fl + 2vr\ 
£ = ( — 3 — ) ^ ^ 

£ = + A sin I — - — I (A.13) 

with 




A — \ — (A.14) 



9 = arcsin^-^ (A.15) 

The first of these energy expressions evaluates to the middle value of Ei, E 2 , and E 3 
when y = 0. When y = 0, the second evaluates to the maximum of Ei, E 2 , and E 3 , and 
the third expression evaluates to the minimum of the three unperturbed energies. 

Appendix A. 4- Algebraic expressions are inconvenient 

It is possible to develop algebraic expressions by taking advantage of the addition 
formula 

sin(a + b) = sin a cos b + cos a sin b (A. 16) 
The energy eigenvalues then can be written as 

E = El + E2 + E3 +Asin (t) (A.17) 



3 13 



E . . Ei + E 2 + E 3 ^ 



(A.18) 
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(A.19) 



To proceed, we require explicit expressions for sin(0/3) an d cos(6 l /3). We begin by 
writing for sin 6* 



sin6> 



e' ie - e~ 



2i 



Zip 
2a§ 



We can solve for e %e to obtain 



e = 



.3§/3 

1 3 



± i l- 



27/3 2 
4a 3 



(A.20) 



(A.21) 



Either choice of sign is acceptable, but we will choose a + sign for what follows. This 
allows us to write 



sm|- 



cos,- 



1 

2i 

1 

2 



/ 27 '/3 2 32/3 

1 : — s % o- 



4a 3 



2a2 



27/? 2 35/3 

Z q~ 



4a 3 



2a5 



i. \ 

" 3 



27 (3 2 .35/3 



+ 



27/3 2 .32/3 
4^ 3__Z 2^i 



> (A.22) 



> (A.23) 



Unfortunately, these expressions are complicated in a way that makes them inconvenient 
for calculations. 



Appendix A. 5. Analytic WKB expression 

It is possible to use these results to develop analytic expressions for the WKB 
approximation for the energy levels discussed in Section 3. In the case of the middle 
level, we may write the approximation as 

Bhbn,!*) = - f -S^fdy (A.24) 

7T J-e V e — V 

After substituting in for E 2 (y) we obtain 



E 2 (gi,g2) 



Ei + E 2 + £*3 



1 



+ 



sin - arcsin 
1 r V3 



my) 



-dy 



7T J-e y/e — y 2 

Similar analytic expressions can be written for the other two levels directly. 



(A.25) 



Multiphoton Bloch-Siegert shifts and level-splittings in a three-level system 



29 



References 

[I] Bloch F and Siegert A 1940 Phys. Rev. 57 522 
[2] Shirley J 1965 Phys. Rev. 138, B979 

[3] Cohen-Tannoudji C, Dupont-Roc J, and Fabre C 1973 J. Phys. B: At. Mol. Phys. 6 L214 

[4] Pegg D T 1973 J. Phys. B: At. Mol. Phys. 6 246 

[5] Ahmad F and Bullough R K 1974 J. Phys. B: At. Mol. Phys. 7 L275 

[6] Hattori T and Kobayashi T Phys. Rev. A 35 2733 

[7] F0rre M 2004 Phys. Rev. A 70 013406 

[8] Ostrovsky V N and Horsdal-Pedcrsen E 2004 Phys. Rev. A 70 033413 
[9] Graham H and Hohncrbach M 1984 Z. Phys. B 233 
[10] Ciblis M B et al 1991 J. Phys. A: Math. Gen. 24 1661 

[II] Wang K, Ho T and Chu I 1985 J. Phys. B: At. Mol. Phys. 18 4539 
[12] D'Andrea A Phys. Rev. A 39 5143 

[13] Matisov B, Mazets I and Windholz L 1995 Quantum Semiclass. Opt. 7 449 
[14] Liu Z, Lin Y, Shang K and Zeng L 1999 Physics Letters A 264 137 
[15] Klimov A B, Sainz I, and Chumakov S M 2003 Phys. Rev. A 68 063811 
[16] Radmorc P M and Knight P L 1981 J. Phys. B: At. Mol. Phys. 15 561 
[17] Li X, Lin D L and Gong C 1987 Phys. Rev. A 36 5209 
[18] Cardimona D A 1990 Phys. Rev. A 41 5016 
[19] Wu Y and Yang X 1997 Phys. Rev. A 56 2443 

[20] Klimov A B, Sanchez-Soto L L, Navaroo A and Yustas E C 2002 J. Mod. Optic. 49 2211 
[21] Bougouffa S and Kamli A 2004 J. Opt. B. Quantum Semiclass. Opt. 6 S60 
[22] Abdcl-Wahab N H 2007 Phys. Scr. 76 244 

[23] Georgi H 1999 Lie Algebras in Particle Physics (New York:Perseus Books) p 98 
[24] Hagelstein P L and Chaudhary I U, to appear in J. Phys. B; available as Preprint quant- 
ph/0709.1961 

[25] Hagelstein P L and Chaudhary I U, to appear in J. Phys. B; available as Preprint quant- 

ph/0709.3557 
[26] Wagner M 1979, Zeit. fur Physik B 32 225 
[27] Yoo H I and Eberly J H 1985 Physics Reports 118 239 
[28] Larson J and Stcnholm S 2006 Phys. Rev. A 73 033805 



